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Despite many measurements on the bulk proper- 
ties of rotating He II, the first local probing of 
the velocity field occurred only in 1960 when Pel- 
lam’ studied the torque on a Rayleigh disk im- 
mersed in a rotating cylinder filled with He II as 
a function of disk position and of temperature. The 
most surprising result of his study occurred in the 
torque versus temperature measurements, which 
seemed to indicate that only the superfluid com- 
ponent of He II exerted torque upon the system, 
so that just below the lambda transition the torque 
vanished entirely. This result appears to be in 
conflict with all theories of He II, and provided 
the motivation for the present study in which the 
velocity field was studied in a rotating system 
using a different type of detector, which measures 
the drag forces exerted on disks placed normal to 
the stream lines of the rotating helium. 

The detectors consisted of pairs of glass disks 
mounted on quartz crossbeams and suspended in- 
side a rotating cylinder (inside diameter 5.3 cm) 
by means of quartz fibers fixed at their upper ends 
in the laboratory coordinate system. The rotating 
cylinder was contained in a partially silvered glass 
Dewar. The bottom of the rotating cylinder was 
sealed by a brass plate, the only connection to the 
main helium bath being through two 0.005-in. fill- 
ing holes drilled in the plate about 0.5 mm off axis. 
The temperature was determined from the helium 
vapor pressure in the cylinder (scale T55p), there- 
by assuring that no errors were introduced from 
flowing He gas. Drag forces on the disks resulted 
in deflection of the detectors through angles which 
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could be measured by means of a light beam re- 
flected from a mirror mounted on the detector, 
and thence onto a circular translucent scale 28 
cm from the axis of rotation. A range of repro- 
ducible angular rotational velocities from 0.21 to 
6 rpm was produced by a synchronous motor cou- 
pled to speed reduction gear boxes and to a lathe 
gear box. A leather belt decoupled the drive from 
the experimental Dewar, thereby reducing vibra- 
tion to a negligible level. 

Two detector assemblies were used. In assem- 
bly I the detectors were 1.2 cm in diameter and 0.2 
mm thick. They were mounted on a diameter at 
2.5 cm center-to-center separation, and suspended 
from a quartz fiber of torsion constant 8x10~ 
dyne cm/rad. In assembly II the detectors were 
0.6 cm in diameter and 0.2 mm thick, mounted at 
a center-to-center separation of 2.3 cm. The tor- 
sion constant was 4.7107 dyne cm/rad. The fi- 
ber constants were determined from the observed 
deflections in air at 300°K at low Reynolds num- 
bers where Stokes’ law obtains, so that the drag 
forces were linear in velocity.2, The constants ob- 
tained in this way were in reasonable agreement 
with those found by timing the fibers as torsion 
pendula using a known moment of inertia. An ad- 
ditional advantage of the air calibration procedure 
is that errors arising from the finite size of the 
disks tend to cancel when deflections in air are 
compared to those in helium. System I was addi- 
tionally calibrated in CS, and in liquid N,. System 
II was calibrated in He gas at 75°K and 580 mm Hg 
and at 2°K and 23.7 mm Hg as well as in air. The 











331 







VOLUME 7, NUMBER 9 







PHYSICAL REVIEW LETTERS 


NOVEMBER 1, 196] 








30 T T 


a | 
STOP ROTATION 
26+ T= 2.153°K - 4 














E 26+ | v=1.5 RPM a 
-) | 
Sag) | 12 em DISKS J 
3 j 
Bol. | RotaTion 
ROTA 

20} | ¥ 
w 18+ 4 
Oo 

16 ’ 

14 

12 

0 12 24 36 48 60 





t (MINUTES) 


FIG. 1. Development of drag force on a disk in 
rotating He II. Rotation is begun at 4 minutes and 
stopped at 36 minutes. 


calibrations were self-consistent to within about 
+20 %. 

Three types of measurement were made with 
these detectors: 

(1) For both systems the buildup of the drag 
forces with time after rotation was begun followed 
the same type of curve and had the same time 
scale as was found in the Rayleigh disk experi- 
ments.' That is, when rotation was started from 
rest a time delay was observed before any drag 
was seen, and the equilibrium drag was attained 
only after about $ hr. When rotation was stopped, 
the drag force immediately decreased, and the 
zero position was regained in a time much shorter 
than the buildup time. A typical drag versus time 
curve is shown in Fig. 1, for system I at a tem- 
perature of 2.153°K and 1.5 rpm. Rotation was 
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begun at ¢=4 minutes and stopped at ¢=36 minutes, 
Thus in this experiment as well as in the Rayleigh 
disk experiments, there is evidence for a velocity 
front propagating from the walls upon initiation of 
rotation, and there is rigidity associated with the 
liquid upon cessation of rotation. 

(2) The temperature dependence of the equili- 
brium drag force was studied in both systems at 
fixed rotational velocities. The results for system 
I are shown in Fig. 2 for a speed of 4.8 rpm. The 
drag force is seen to be almost independent of 
temperature from 1.3°K to the lambda tempera- 
ture, 7), with a slight minimum near 1.9°K. The 
drag appears to be continuous across the lambda 
transition. However, above T) fluctuations in the 
detector became quite large, so that it was very 
difficult to take measurements and the points are 
subject to large error. The constancy of the drag 
force with temperature is in striking contrast with 
the torque measurements obtained with the Ray- 
leigh disk. The temperature variation of drag 
obtained with system II at fixed rotational veloci- 
ties is similar to that for system I. 

(3) At fixed temperatures below T) the drag 
force was studied as a function of the velocity (v) 
of rotation. For system I the drag force was 
found to be quadratic in velocity at all tempera- 
tures not too near T,. This is the velocity de- 
pendence expected for a disk immersed in a 
classical fluid at Reynolds numbers & greater 
than about 50.” In the present experiments, 

R (=pvd/ny,) always exceeds 100 so that classical 
drag may be expected. The dimensionless drag 
coefficient Cp (defined as the ratio of the observed 
drag force to pv? times the area of the disk) ob- 
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FIG, 2. Temperature dependence of the drag force on a disk in rotating helium. 
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tained from the data of Fig. 2 is about 0.25 at all 
temperatures. This drag coefficient is approxi- 
mately equal to that obtained in the sphere drag 
measurements of Dowley, Firth, and Hollis Hal- 
lett,s an agreement which is somewhat surpris- 
ing since the classical drag coefficient in this 
range of Reynolds numbers is 0.47 for spheres, 
while it is 1.17 for disks.? The behavior in the 
smaller system II was much different from sys- 
tem I. Except close to T) the drag force at fixed 
temperature varied as v” with 1 <n<1.3. Figure 
3 presents some representative drag versus rota- 
tional velocity measurements for systems I and 
I. At temperatures within 0.1°K of T) the be- 
havior of both systems becomes more compli- 
cated. The exponent m becomes a function of ve- 
locity, but always approaches two at high veloci- 
ties. The high-velocity drag coefficient is always 
about 0.25. 

The significance of the linearity of the drag 
force with velocity in system II is not yet clear, 
and may well be an instrumental effect. This be- 
havior and the behavior close to T, are still under 
study. However, because of the present lack of ex- 
planation for the Rayleigh disk results, we venture 
the following suggestion which offers a consistent 
method for interrelating those results with the 
linear drag found in the present experiment. 

In classical fluids drag forces linear in velocity 
occur only for flow at low Reynolds numbers. The 
linearity of the system II measurements and the 
instantaneous decrease in the drag force upon 
cessation of rotation suggest that there is rigidity 
associated with rotating He II. This rigidity could 
arise from the vortex lines which must exist in 
rotating He II, and which probably exert drag 
forces upon bodies placed in rotating He II. Since 
the vortex lines are excitations of ground-state 
superfluid, they must be associated with the nor- 
mal fluid, and would be seen in a drag experiment 
as an excess normal fluid viscosity. The Reynolds 
number describing the normal fluid flow would 
then be quite small, so that Stokes’ drag linear 
in velocity could occur with attendant flow separa- 
tion. Since the vortex line density is independent 
of temperature, the drag force is also. The line 
density is proportional to the rotational velocity, 
but the excess viscosity is not necessarily propor- 
tional to the line density. Such an excess normal 
fluid viscosity should influence the Rayleigh disk 
experiment. The equation for the torque on a 
Rayleigh disk is derived under the assumption of 
high Reynolds number flow with no flow separa- 
tion. In rotating He II the total torque on a disk 
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FIG. 3. Variation of the drag force on disks in Hell 


with rotational velocity. Two distinct types of behavior 
are found. 


with this assumption should be proportional to 
Psv," +P,V,", and presumably v, =v, since rotat- 
ing He II is known to possess classical angular 
momentum.‘ If flow separation of the normal fluid 
occurs, the second term vanishes, and the total 
torque is proportional to p 5? AA in agreement with 
Pellam’s observations. 

Vinen® has argued that vortex line turbulence 
may arise from heat currents in He II. If the 
above speculations are valid, it should be possi- 
ble to increase the excess viscosity by heat cur- 
rents. Preliminary experiments tend to support 
this idea, but large fluctuations of the detector 
due to unstable flow prevent definite conclusions 
at this time. 

The author would like to thank Professor J. R. 
Pellam and Professor J. G. Dash for valuable 
discussions, and Mr. Donald Clinton and Mr. Earl 
McGowan for the construction of the apparatus, 
and for continuing advice on its design. 
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Cohen and Falicov' have discussed the effect of 
spin-orbit coupling on the topology of the Fermi 
surfaces of the hexagonal close-packed metals. 
They conclude that spin-orbit splitting lifts the 
degeneracy between the first and second bands at 
the corners of the Brillouin zone and separates 
the small pockets of holes (arm caps) in the first 
band from the region of holes (arms) in the second. 
This gives rise to a closed-hole surface (cap sur- 
face) in the first band and a second-band surface 
that is infinitely extended along the hexagonal axis 
[ Fig. 2(b) of reference 1] due to the intersections 
of arms. Both of these surfaces have extremal 
cross sections at the hexagonal face of the zone 
which one would expect to be observable in de 
Haas~—van Alphen, cyclotron resonance, and ul- 
trasonic attenuation experiments. 

During the course of detailed de Haas—van Al- 
phen studies of zinc and cadmium, we have ob- 
served oscillations in the magnetic susceptibility 
whose periods correspond to estimated values for 
the aforementioned extremal cross sections and 
which we have been unable to account for with 
Fermi surface models that neglect spin-orbit 
coupling. The de Haas—van Alphen periods were 
obtained by a torsion method in magnetic fields 
up to 23 kilogauss, using a null deflection tech- 
nique and automatic recording. The angular de- 
pendences of the periods essential to this dis- 
cussion are shown in Fig. 1. @ is the angle be- 
tween the magnetic field direction and the hexag- 
onal axis. 

In Cd, for small @ the values of the A and x 
periods can be determined to within approximately 
2% because all of the oscillations present are of 
comparable amplitude. As @ is increased the 
amplitude of the main arm oscillation increases 
rapidly, while effects due to arms in the other 
symmetry directions [periods indicated by dashed 
lines in Fig. 1(b)] are hardly observable. The 
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amplitudes of the x and A oscillations drop rapidly, 
with A undetectable at 6=1° and x at @=10°. 

In Zn, analysis is more difficult because x and 
4 occur only as a weak modulation on a strong 
carrier—the arm periods. It was possible to ob- 
tain an unambiguous value for the A points out to 
a 6 of 6°. However, there is no way of distinguish- 
ing from our data whether x lies below the arm 
period as shown in Fig. 1(a) or an equivalent dis- 
tance above. The lower value was selected to be 
consistent with the relative sizes of the spin-orbit 
splitting in Zn and Cd (reference 1) and with the 
effect of this splitting on the difference in cross 
section of arm junction and cap surface as ex- 
plained below. For @ greater than 7°, beating 
among the three arm periods makes detection of 
the A oscillation difficult due to its decreasing 
amplitude and to the presence of a comparable 
arm period. Therefore, the angular extent of the 
A period is not sharply determined. Near @=15, 
where the A period might beat distinctly with the 
arm periods, a careful search yielded no sign of 
its presence up to 23 kilogauss nor did any unex- 
plained periods appear at larger 6. 

The arm data in Zn agree with the few periods 
observed by Verkin and Dmitrenko,? while those 
in Cd fit very well with the work of Berlincourt,; 
but no mention of the x and A periods appears in 
either paper. We feel that the identification of 4’s 
as cap surface periods and x’s as arm junction 
periods is correct because the other periods ob- 
served in Zn and Cd agree with the general form 
of portions of the Fermi surface predicted by the 
single orthogonalized plane wave construction,‘ 
and because this construction requires arm junc- 
tions and cap surfaces in the single zone scheme. 

The presence of the spin-orbit energy gap de- 
creases the cross section of the cap surface and 
increases that of the arm junction from the com- 
mon value that these cross sections would have in 
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FIG, 1. Angular dependence of the arm periods (0) 
and the periods attributed to the arm junctions (x) in 
the second band and the cap surface (A) in the first 
band. (a) Zinc with the magnetic field lying in a (1120) 


plane (in reciprocal space). (b) Cadmium with the magnetic field lying in a (1010) plane (in reciprocal space). The 
dashed lines are drawn for the other arm periods expected using cross sections described by the full line. 


the limit of a vanishing gap. If we use the nearly 
free electron approximation to estimate the in- 
fluence of the size of the energy gap on the dif- 
ference between cap surface and arm junction at 
the Brillouin zone face, the periods of Fig. 1 for 
magnetic field along the hexagonal axis yield a 
ratio of the spin-orbit energy gap in Cd to that in 
Zn of 2.2. This value is in reasonably good agree- 
ment with the estimate of 3 given in reference 1. 
The alternative possible assignment for the x 
period in Zn, discussed above, leads to an energy 
gap in Zn greater than that in Cd, an unlikely situ- 
ation. 

Our observation of periods due to the cap surface 
and arm junctions in Zn in magnetic fields as high 
as 23 kilogauss seems at first glance surprising in 
view of the possibility, first pointed out by Cohen 
and Falicov,® of “magnetic breakdown” of the spin- 
orbit energy gap. Blount® has shown that the rele- 


vant parameter determining breakdown is hw,E f/ 
E g3 i.e., when this quantity is much greater than 
one, the semiclassical electron orbits are deter- 
mined by ignoring the energy gap in question. 
(Here w, is the cyclotron frequency, Eg is the en- 
ergy gap, and Ep is the Fermi energy.) Using 
Cohen and Falicov’s value for E g in Zn,' we obtain 
for a field of 23 kilogauss a value for the above 
parameter that is of the order of one. Thus one 
would expect behavior characteristic of a transi- 
tion region with additional scattering processes 
which would seriously decrease the amplitude of 
de Haas—van Alphen oscillations. 

However, the criterion for breakdown should 
depend upon the inclination of the semiclassical 
orbit relative to the Brillouin zone plane under 
consideration. Intuitively, it would seem that if 
the electron’s orbit cuts the Brillouin zone plane 
at a steep angle, then the probability for “jumping 
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the gap” ought to be higher than if the electron 
moves parallel to the zone plane. This result is 
borne out by detailed calculations; in fact, when 
the electron moves parallel to the zone plane, 
magnetic breakdown is completely inhibited. Thus, 
the spin-orbit energy gap across the hexagonal 
zone face in such metals should be seen, even at 
high fields, when the magnetic field is directed 
along the hexagonal axis of the crystal. 

We have calculated the transition probability for 
scattering out of a semiclassical orbit as a re- 
sult of the perturbation: V=V,cos(K,r). Accord- 
ing to nearly free electron theory, E gm 2V9. The 
method used is an alternative one to Blount’s and 
is one suggested by Cohen.” We start with the 
transition probability per unit time due to scatter- 
ing: 


w=(20/AX>, ,1V,,,06,,-E,)),, 


where ( ), indicates an average over the orbit out 
of which the electron is being scattered. When the 
plane of the electron’s orbit makes a small dihe- 
dral angle @ with the Brillouin zone plane, we ob- 
tain for the transition probability per cyclotron 
period: 


7 = 2 
w/w, ‘. /32 hw OE 


rp sin’o, (2) 
where ¢ is the polar angle of the electron’s orbit 
(measured with respect to the magnetic field). 
When the orbit cuts the Brillouin zone plane at 
right angles, the calculation of w/w, may also be 
effected easily; the result for this case differs 
from (2) primarily through the deletion of the 
factor 6. By a WKB argument Blount® has shown 
that when w/w, >1, the probability for electron 
tunneling is exp(-7w/8w,). Combining Eq. (2) with 
Blount’s exponential factor and the E g of reference 


1 suggests that the effects of spin-orbit coupling 
should be observable within an angular range of 
about 6° from the hexagonal axis in zinc. 

One would expect the arm junction period to 
exist over only a limited angular region due to the 
disappearance of an extremal cross section for 
this region of the Fermi surface, but the cap sur- 
face period should be present (barring magnetic 
breakdown) over the entire range of 6. We have 
observed the cap surface period for 6 as large as 
6° in Zn and 1° in Cd, and the arm junction period 
for 6 as large as 1° in Zn and 10° in Cd. It is ob- 
vious from this that our data to date can cast little 
light on the question of the quantitative correct- 
ness of the present theory of magnetic breakdown, 

We have been unable so far to observe cap sur- 
face and arm junction periods in magnesium. 

We wish to thank J. W. Rutter and B. W. Roberts 
of the General Electric Research Laboratory at 
Schenectady for generously supplying the cadmium 
single crystal used in this work. 
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In a recent experiment,’ the magnetic flux 
trapped within a doubly connected superconductor 
has been measured by cooling a thin ring embrac- 
ing an external flux below the transition tempera- 
ture. The experimental result is plotted in Fig. 1, 
which gives the total flux ¢ through the ring (which 
can never be changed in the superconducting state) 
as a function of the external flux ¢,. The quantity 
E also plotted in Fig. 1 will be discussed later on. 

The unit ¢,=hc/e of Fig. 1 is the unit of quanti- 
zation briefly mentioned by London? and by Onsag- 
er? which stimulated the experiment. In order to 
find an indication for the quantization in units of 
o,/2, let us first discuss the argument of London. 
After this we will give a simple argument that the 
experiment gives an indication of the existence of 
bound pairs in the superconductor. 

If Y,(j"++Ty) is the wave function of the ground 
state of the superconductor, the wave function 


y= expliDS&,)/nly (1) 
corresponds to a superfluid flow with a velocity 
field v(t) given by 

. DF) = grad S(F)=mv(F)+ /c)A@), (2) 
thus obeying London’s equation, 
curl[mv¥+ (e/c)A]=0. (3) 


For a closed curve C located entirely within the 
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FIG. 1. The experimental values of ¢ (step-like 


curve, schematically) and the kinetic energy E of the 
Superconductor as function of ¢¢. 
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superconductor and not embracing any hole, the 
integral 


$B-d8=$_,[mv+ (e/c)A]-d8=(e/c)> (A) 








therefore must be zero. This, however, need not 
be the case if C embraces a hole. As a conse- 
quence of (3) the integral then has the same value 
for all curves around the hole. If the ring is thick 
enough, the curve can be placed in its interior so 
that the current density ~v(r) can be neglected. In 
this case one has = /Hdf. 

The quantization of the flux ¢ then is supposed 
to follow from the fact that the wave function is 
single valued. London states, without proof, that 
this requires $p-ds=nh. This sounds plausible 
since y, is single valued so that exp[i>, S(r;)/h] 
also has to be single valued in order to guarantee 
the same for y. 

We now want to show that this is not necessary 
if y, contains strong correlations of the kind pro- 
posed by Bardeen, Cooper, and Schrieffer.* 

According to BCS the ground-state wave function 
can be written as 



















where A indicates the antisymmetrization proced- 
ure and w(f) is the Fourier transform of v,/ux 
=w} (ux, and v;, being defined as in BCS). 

Let us now introduce cylinder coordinates with 
the axis along the axis of the ring. If we do not 
indicate the y and z dependence explicitly, (5) re- 
duces to 







= Allw(g, - 9;). (6) 





¥9 
As a first approximation one may expect 

w(y, - ¢,) to be similar to the wave function in 
an infinitely extended homogeneous superconduc - 
tor, i.e., only different from zero in the neighbor- 
hood of y,=¢,, which would be single valued in 
(y, + @)/2 but not in g, and gy, individually. The 
latter would require 


w(¢y)=w(g+ 2m). (7) 


This condition, however, can be fulfilled in re- 
placing w by 












*” 2 wy +27v), (8) 







One also expects excited states corresponding 
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to pairs running around the axis with a total angu- 
lar momentum n of the type 
in (9, + P2)/2 * 

e w Go ?,)- (9) 
Here n has to be an integer in order to make (9) 
single valued in (y,+¢,)/2. The single valuedness 
in g, and ¢, individually can be saved if Wy is re~ 
placed by 


w =D, w(p+2nv)\(-1)", (10) 


if m is odd. 

The wave functions (9) and (10) are not exact 
eigensolutions of the Schrédinger equation, but 
they are good approximations deviating from the 
solutions of the homogeneous medium in only a 
very small part of the configuration space 
0 <¢,, ~,< 27 if the extension of the pairs is small 
compared to the dimensions of the ring. The wave 
functions of the whole system are given by 


ind). ,/2 
e Allw .- g.) for n even, 
2% %;) 


ind).¢./2 
«| lUFR w 9, - ?;) for n odd. (11) 


We do not consider excitations corresponding to 
the breakup of pairs and the acceleration of a 

pair relative to the rest of the pairs since these 
are only excited by fields above the critical field. 
The application of (4) then gives the required half- 
integral values of e@/hc. 

The discussions above only give the wave func- 
tions which are compatible with the requirement 
of single valuedness but do not tell us which of 
them actually is obtained for a given external flux. 
For this purpose let us consider the energy of the 
electrons as a function of the external flux ¢,. 
Since there is no first order change of the wave 
function of a superconductor in an external field, 
the energy of the excited states above the ground 
state is just the kinetic energy of the superfluid 
flow which can be estimated by assuming v(r) to 
be constant inside the penetration depth and zero 
outside. This gives as the excitation energy per 
particle 





E/N =tmv*(,/Ar), (12) 


where Ar is the thickness of the ring and \ the 
penetration depth. 
The velocity v then can be calculated from Max- 


well’s equation: curl H| = 4in,ev/c (n, is the num- 


ber of superconducting electrons). If one replaces 
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1v*H by the flux (¢ - og) produced by the supercon. 
ducting electrons, one obtains finally 


n.ev=(o- o,)e/407n; E =const(¢ - >). (13) 


The result is plotted in Fig. 1 for different values 
of the total flux ¢. 

Whereas in the superconducting state the transi- 
tions between different parabolas, even if energet- 
ically possible, are forbidden because of d¢/dt=0, 
one expects that the state obtained by cooling the 
ring down from the normal state with an external 
flux will be the state of lowest possible energy, 
i.e., that state which can be obtained by the small- 
est possible acceleration of the electrons. This 
would explain completely the experimental result 
of Fig. 1, in particular the fact that the first flux 
quantum is trapped already if the external flux has 
only reached the value of ¢,/4. 

The model considered above is of course very 
poor. We have only treated the case of zero tem- 
perature and we have left out all the complicated 
events which happen during the transition from the 
normal to the superconducting state. Nevertheless 
we hope that the indications are strong enough to 
believe that the experiment of Doll and Nabauer is 
a macroscopic) consequence of angular momentum 
quantization (similar to the quantized vortex lines 
in He II) and is directly related to the existence of 
bound pairs of electrons in the superconductor. 

Note added in proof. Meanwhile there has ap- 
peared a Letter by Onsager® explaining the occur- 
rence of odd multiples of hc/2e in essentially the 
same way as above. We can use our wave function 
to exhibit the equivalence of this idea to the pairing 
condition as used by Byers and Yang.® For this 
purpose one only has to notice that introducing 
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The wave functions (11) then can easily be trans- 
formed to the usual notation in second quantization: 
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Inside a thick ring the kinetic angular momenta 
are given by h(v - e¢/hc) if hv is the canonical 
angular momentum. Therefore with ¢=hcm/e 
(o=hc(m +4)/e] the pairing with zero current 
occurs for (v, -v)[v+4, -v-4)] in the case of 
even (odd) n. The description in terms of wave 
functions (11) or (15) has the advantage of not 
only showing the equivalence of two apparently 
different pictures but also of being valid in the 
case of thin rings where the current does not van- 
ish as assumed by Byers and Yang. 

We wish to thank H. P. Duerr, W. Heisenberg, 
H. Koppe, M. Nabauer, A. de-Shalit, and V. F. 


Weisskopf for stimulating discussions. 
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OSCILLATORY GALVANOMAGNETIC EFFECT IN METALLIC SODIUM" 


R. Bowers, C. Legendy, and F. Rose 


Laboratory of Atomic and Solid-State Physics and Department of Physics, Cornell University, Ithaca, New York 
(Received October 5, 1961) 


We wish to report the observation at 4°K of a 
novel galvanomagnetic effect in sodium of very 
high conductivity (p399°K/p4°K = 7500; w.7 = 40 
in 10000 gauss). In this effect, the Lorentz force 
produces an electrical oscillation of low frequency 
(ranging from 5 to 35 cycles/sec in our experi- 
ment). Our results provide direct evidence of a 
particular macroscopic electron gas excitation 
described by Aigrain' and given by him the name 
“helicon” excitation. There is also a connection 
between our observation and certain effects ob- 
served in bismuth at microwave frequencies.” 

The experimental arrangement (Fig. 1) is based 
on a configuration for measuring resistivity due 
to Bean, DeBlois, and Nesbitt.* A cylindrical 
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FIG, 1. Schematic diagram of experiment. Primary 
coil ~300 turns of No. 28 copper wire; secondary coil 
*1200 turns of No. 43 copper wire. Specimen and coils 
immersed in liquid helium at 4.2°K. 


specimen of sodium (height 30 mm and diameter 
4 mm) is placed inside primary and secondary 
coils. On closing or opening the primary circuit, 
eddy currents are induced in the specimen which 
can be detected by the secondary coil. From the 
time constant of the observed voltage in the sec- 
ondary, one can estimate the resistivity of the 
specimen. For a sufficiently long time after the 
primary pulse, the decay is a simple exponential. 
Our experiment differs from that in reference 3 
in that a constant magnetic field (H,) is applied 
perpendicular to the axis of the specimen and 
coils. The constant field H, is much larger than 
the time-dependent magnetic field (h = few hundred 
gauss) produced by the primary coil. The voltage 
of the secondary coil is observed on closing or 
opening the primary circuit. Figure 2 shows 
oscilloscope traces of the secondary voltage as 
a function of time for various values of the bias 
field H,, using sodium for which p3q9°%/p4°K 
~ 7500 (mean free path at 4°K=4 mm). In zero 
magnetic field the normal eddy current decay is 
observed. The other three traces show the effect 
of applying H,. One sees damped oscillations 
whose frequency is proportional to magnetic field. 
The frequency at 10000 gauss is 32 cycles/sec. 
We believe the oscillations to originate as fol- 
lows. The abrupt change in the primary field (2) 
excites stable modes of the electron gas in which 
the electric vector is circularly polarized about 
the large magnetic field.’ The oscillation observed 
is the component of this circular wave resolved in 
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FIG, 2. Voltage in secondary coil as a function of 
time on opening primary circuit for different values of 
Hy. Starting from the top, the four curves correspond 
to Hy=0, 3600, 7200, and 10800 gauss, respectively. 
Ordinate is 5 millivolts per large division; abscissa is 
50 milliseconds per large division. Specimen proper- 
ties: P3o9°K/P4°K ~7500; wT ~ 40 in 10000 gauss. 


the direction of the axis of the secondary coil. A 
precise theory of this effect requires the solution 
of the field equations with the complicated bound- 
ary conditions used in our experiment. However, 
a useful physical picture of the origin of the circu- 
larly polarized waves can be derived from the 
following simplified theory for the infinite solid 

of very high conductivity. 

In a strong magnetic field (w.7> 1), the current 
vector and the electric field vector are nearly per- 
pendicular, as the component of the electric field 
along the current direction becomes very small 
compared to the transverse (Hall) field. We there- 
fore write* 


E+p 0” D x H= 0, 

T=nev pe (1) 
using rationalized mks units and conventional no- 
tation. We use Eq. (1) for fields and currents 
varying in space and time, assuming that the mean 
free path and the relaxation times are short com- 
pared to the corresponding periods in the solutions. 
From Eq. (1) and Maxwell’s equations (neglecting 
displacement current) it follows that 


8E/8t= (1/ne)[V x(V xE)]xH,, (2) 
assuming that H= Hy +h(r, t), iH! > (hi, and H, is 
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along the z axis. We solve by substituting 

E= E, exp|i(w7 - K-r)]. Equation (2) reduces toa 
set of three linear homogeneous equations for the 
components of E,. The secular equation yields 
two _— frequencies for any given wave vec- 
tor k: 


w=tkk H)/ne. (3) 


The corresponding solutions, omitting a common 
multiplicative constant, are: 


7 o« 2 2 ; 
EY, eg. +k, bp em, 0), 


h.= ; 2 2 . 
Ay &. /w U*k, + ik - r. +k, : “2. + ikk ), 
i, = Fkh,. 


For propagation along the z direction, these re- 
duce to k#(1,+i,0); (k°/wu,)(Fi, 1,0); #kh,, respec- 
tively. All the above solutions describe circularly 
polarized plane waves, the two signs correspond- 
ing to right-hand and left-hand screw spatial con- 
figurations. The vectors rotate about the magnetic 
field lines always in the direction corresponding 
to the cyclotron rotation of the carriers. The E 
field in general has a component along the propa- 
gation direction; the h and i fields do not. 

At 10000 gauss the observed frequency was 32 
cycles/sec; with these values Eq. (3) yields a 
half-wavelength of 3.4 mm for a wave that propa- 
gates along the field lines. This is comparable 
with the corresponding specimen dimension which 
is 4mm. The observed linear dependence of the 
frequency on the magnetic field is also in agree- 
ment with Eq. (3). 

We have not been able to find any other explana- 
tion for the observed oscillations. Mechanical 
oscillations of the coil in the magnetic field are 
expected to produce an oscillation whose fre- 
quency is determined by the mechanical spring 
constants and is independent of H,; the amplitude 
of the resulting oscillations would depend on H, 
in a manner quite different from that shown in 
Fig. 2. Moreover, we have taken care to reduce 
mechanical effects to the minimum. The oscil- 
lations observed occur at frequencies which are 
too low by orders of magnitude to be associated 
with Landau levels. 

The simple theory predicts that the frequency 
should be independent of conductivity in the high- 
conductivity limit. This was verified on investi- 
gating the oscillation in four sodium specimens 
whose residual resistance ratios were 7500, 4500, 
3200, and 1200, respectively. It was found that 
oscillations were more damped as the conductivity 
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decreased but the frequency did not change appreci- 
ably. At a ratio of 1200, no oscillations were 

found because of excessive damping, indicating 

that the conductivity must be very high in order 

to see the oscillations. Sodium is very favorable 

in this respect since it has a very small magneto- 
resistance. 

In a variation of the experiment described above, 
the oscillation was excited under resonant condi- 
tions. Helmholtz coils were mounted (outside the 
cryostat) so that their axis was perpendicular both 
to the large field H, and to the axis of the second- 
ary coil. The latter orientation reduces direct 
mutual inductance coupling. In Fig. 3 is shown the 
response of the secondary to a sinusoidal excitation 
of the Helmholtz coils; the specimen was the same 
as that used for Fig. 2. The amplitude of excitation 
current was maintained constant and the frequency 
varied. A resonant response was obtained. From 
studies of different specimens, it was found that 
the Q for the resonance increases with increasing 
conductivity but the resonant frequency is nearly 
independent of the conductivity in the range investi- 
gated. As yet, no identifiable harmonics have been 
found. Although the dependence of the resonance 
on specimen size has not been investigated in detail, 
it is observed that the resonant frequency increases 
as the length along H, is decreased. 

We expect to carry out experiments in which the 
geometrical arrangement is more amenable to de- 
tailed theoretical analysis than the present arrange- 
ment; this requires modification of our magnet- 
cryostat system. This work provides further evi- 
dence that one can study in a solid some of the 
macroscopic electron gas excitations more fre- 
quently investigated in gaseous plasmas.? 

We are indebted to J. Krumhansl and M. T. Tay- 
lor for helpful discussions, and to S. Tallman and 


FIG. 3. Voltage in secondary as a function of fre- 
quency of exciting coil for fixed H)=10800 gauss. Or- 
dinate is 1 millivolt per large division; abscissa is 5 
cycles/sec per large division. Specimen was the same 
as that used for Fig. 2. 


E. Johnson for contributions to the experiment. 
The high-purity sodium was prepared by C. E. 

Taylor of the Lawrence Radiation Laboratory, 

Livermore, California. 
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CRITICAL POINTS AND FERROMAGNETISM 


E. P. Wohlfarth and J. F. Cornwell* 


Department of Mathematics, Imperial College, London, England 
(Received October 16, 1961) 


Ordinary and some singular critical points of 
E(k) in the Brillouin zones of metals produce infi- 
nite discontinuities of N’(E), the derivative with 
respect to energy E of the density of states N(E).'»” 
As a result, N(£) curves contain a number of 
sharp and narrow peaks, the structure being more 
pronounced for the more complicated metals. This 
is illustrated in Fig. 1, showing N(E) for body- 
centered iron. The curve was calculated, as de- 
scribed by us previously,*® by applying the Slater - 
Koster interpolation scheme‘ to the results of 
Wood,* who used the augmented plane wave method 
but did not, according to the present calculations,’ 
cover a sufficient number of points in the zone to 
give a fully reliable N(EZ) curve, nor did he inves- 
tigate the critical points. 

The N(E) curve of reference 3 was obtained by 
finite difference methods and thus also did not ex- 
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hibit the sharp peaks of the curve of Fig. 1; this 
was calculated by numerical differentiation. Fur- 
ther consideration of this curve has now become 
of interest, for the following reason: The peaky 
structure, here exemplified for iron but present 
for all metals,®° seems to have a direct bearing on 
the occurrence of ferromagnetism, and may thus 
explain the recently discovered ferromagnetism 
in the Sc-In alloy system’; the properties of iron 
itself may also be discussed on this basis.* 

It is thus suggested that on adding indium to 
scandium the Fermi energy passes through a very 
narrow peak in N(E) corresponding effectively to 
the very narrow composition range 23.8% to 24.24 
In over which ferromagnetic behavior was found to 
occur so sensitively. The low value of the Curie 
temperature T,., observed to be about 6°K for 
these alloys, is then related to the low effective 





FIG, 1. Density-of-states 
curve for body-centered iron. 
N(E), density of states per atom; 
E, energy in rydbergs. 
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degeneracy temperature proportional to the width 
of the peak. The susceptibility x above T e could 
only be calculated if the precise structure of the 
peak and the exchange energy were known. How- 
ever, as exemplified in Fig. 7 of Elcock et al.,° 
and as frequently discussed for normal ferromag- 
netic metals and alloys,° the 1/x, T curves calcu- 
lable on this basis frequently appear to be very 
closely linear over relatively wide ranges of tem- 
perature above 7. The slopes of such linear 
curves then bear no relation whatever to the sat- 
uration magnetizations obtained from low-tempera- 
ture measurements. For the Sc-In alloys there 
was indeed no such relation: The saturation mag- 
netic moments per atom are given as 0.035 to 

0.039 Bohr magneton while the values derived 

from the apparent Curie-Weiss constants ranged 
from 0.65 to 0.70 over the ferromagnetic composi- 
tion range. This divergence shows that a localized 
model is not applicable here, while the discussion 
just given shows that the observations can be ex- 
plained in terms of a band model including a peaky 
structure. Further evidence bearing on this inter- 
pretation would come from low-temperature spe- 


cific heat measurements. 
We wish to thank B. T. Matthias and A. M. Clog- 
ston for helpful discussions of their results. 
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NUCLEAR SPIN STATES IN MOLECULAR CRYSTALS” 


L. B. Borst, S. L. Borst, L. Koysooko,f H. Patel, and E. Stusnick 
Department of Physics, New York University, University Heights, New York, New York 
(Received August 21, 1961; revised manuscript received October 9, 1961) 


Nuclear spin states are commonly identified in 
gases, but are not generally expected to exist in 
condensed phases. Whereas ortho and para forms 
are acknowledged for hydrogen both in liquid and 
solid,’ restrictions in other cases imposed by 
neighboring molecules presumably prevent quan- 
tized rotation with a clear distinction between odd 
and even rotational quantum numbers with asso- 
ciated parallel and antiparallel nuclear spin align- 
ment. 

In the case of molecular hydrogen, the elastic 
neutron scattering cross section for the ortho 
State is large compared to the para state for both 
liquid and gas (30:1 in the case of the gas). Neu- 
tron scattering measurements therefore permit 
the identification of ortho and para states. 

The experiment consisted of measuring the at- 
tenuation of a neutron beam by the sample under 
investigation which was immersed in liquid heli- 
um. Temperatures below 4°K were achieved by 
evacuating the liquid helium and reducing the 
boiling point. The beam intensity as measured 


by a thin BF, counter was observed as a function 

of time with sudden changes in sample tempera- 
ture, holding other experimental parameters con- 
stant. The neutron beam was derived from a re- 
entrant hole in an aqueous moderator immersed 

in liquid helium. A velocity spectrum of the emer- 
gent beam was best fitted by a Maxwell distribution 
having a temperature of 40°K. 

Solid acetylene (HCCH), a linear molecule, was 
cooled, first to 4°K and then to 2°K, and the total 
cross section was followed as a function of time 
(Fig. 1). The cross section changed from an ini- 
tial value of 140 barns to 10 barns per molecule. 
The temperature was rapidly raised to 77°K, after 
which the cross section returned to its original 
value. The time for 50% conversion from the low- 
temperature form (mostly para) to high-tempera- 
ture form (3 ortho-1 para) was about 10 minutes. 

Heat capacity data for water were once inter- 
preted in terms of ortho and para states.* The 
Pauling model* so completely explained the resid- 
ual entropy that the alternate interpretation was 
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FIG. 1. Time-dependent cross section of acetylene. 


forgotten. Comparison of the scattering cross 
section of molecular hydrogen with water for the 
high-temperature mixture (3 ortho-1 para) shows 
comparable cross sections and energy dependence. 
Even though hydrogen-bonded, water would seem 
to have ortho and para modifications.®»* Trans- 
mission experiments similar to those on acety- 
lene confirm this prediction. 

Water has been studied as a neutron moderator 
at cryogenic temperatures.*:? Allowed transitions 
between ortho and para states are possible by 
means of spin exchange scattering. This inelastic 
process constitutes an additional moderation mech- 
anism.® A kilogram sample was cooled to 4°K. 
The neutron spectrum within the moderator, as 
measured by a foil activation technique, indicated 
the existence of a substantial fraction of the neu- 
trons in a Maxwell distribution near the thermo- 
dynamic temperature.® A beam of neutrons was 
removed under geometric conditions which as- 
sured minimum perturbation of the velocity dis- 
tribution function. 

An absorption curve, using 1/v absorber and a 
thin 1/v detector, showed that 20% of the neutrons 
in the beam from a 4° moderator had an effective 
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cross section in indium of 1500 barns. This soft 
component was not present at moderator temper- 
atures of 77 or 300°K. These data are considered 
evidence that 20% of the neutrons reaching 300°K 
enter a thermal distribution near 4°K. All the 
foregoing experiments were conducted immediately 
after cooling the moderator and were conducted 
within a period of less than an hour. 

A velocity spectrum did not, however, clearly 
delineate this low-temperature thermal distribu- 
tion. Plotting the intensity per channel as a func- 
tion of time disclosed an exponential decay with a 
half-life of about 100 minutes (Fig. 2). Under the 
experimental conditions of constant moderator 
temperature and careful monitoring, the intensity 
change could only be identified with a change in 
the scattering power of the moderator. A reduc- 
tion in scattering power would increase the leak- 
age of neutrons from the refrigerated moderator, 
thereby reducing the thermal source intensity. At 
the end of the refrigeration period the thermal 
neutron distribution returned to its original inten- 
sity. A change in intensity by a factor of ten was 
observed in the course of eight hours for neutrons 
between 500 and 1500 psec/m reciprocal velocity. 
The channel covering 0-500 psec/m shows a less 
pronounced time dependence because it includes 
fast and epithermal neutrons. 

These changes in spectral intensity are consid- 
ered evidence of substantial ortho to para trans- 
formation with associated loss of scattering power. 
After para transformation, the neutron flux is so 
low that it is not possible to distinguish a Maxwell 
distribution near the thermodynamic temperature 
from background effects. During the first part of 
the experiment, when ortho molecules are present, 
and earlier measurements showed evidence of a 
Maxwell distribution near the thermodynamic tem- 
perature, velocity selector statistics are inade- 
quate again to demonstrate its existence. 

Water is therefore an ineffective neutron mod- 
erator at 4°K. In the para state no levels exist 
near the ground state so that there is no effective 
mechanism for inelastic scattering and neutron 
moderation. While the ortho state has an appro- 
priate transition for neutron energy loss, the 
lifetime of the ortho molecule for spontaneous 
transformation to para is short. If the tempera- 
ture is such that a significant fraction of the mole- 
cules exist in the ortho state, then its neutron 
moderating properties should be good. 

Nuclear spin states found in water and acetylene 
are expected in many molecular crystals, the ro- 
tational energy levels being shifted to higher tem- | 













FIC 
tral | 
mode 


per: 
teri 
dihy 
ins 
Cry 
mus 
vati 


*R 
and 
th 
Thai 





961 


oft 
or= 
red 
°K 


ity 


n= 
aS 

ns 
ty. 


ent, 


re 











VoLuME 7, NUMBER 9 


PHYSICAL REVIEW LETTERS 





NoveEMBER 1, 1961 








FIG. 2. Time-dependent spec- 
tral distribution from aqueous 
moderator at 4°K. 
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peratures by the crystal potential. Neutron scat- 
tering is an incisive tool at least in the case of 
dihydrogen molecules because of the large change 
in scattering power between ortho and para states. 
Crystal structure determinations using neutrons 
must now be reviewed in the light of these obser- 
vations, particularly in the case of water. 
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NEW DETERMINATION OF THE K*-DECAY BRANCHING RATIOS* 





Byron P. Roe and Daniel Sinclair 
University of Michigan, Ann Arbor, Michigan 


and 


John L. Brown, Donald A. Glaser,f John A. Kadyk, and George H. TrillingT 


Lawrence Radiation Laboratory, University of California, Berkeley, California 
(Received September 28, 1961) 


Present knowledge of the branching ratios 
among the various decay modes of K* mesons 
has come exclusively from emulsion measure- 
ments. The results of the most precise emul- 
sion determinations are presented in the first 
three columns of Table I. Since some of the ex- 
perimental values of the branching ratios are in 
disagreement by substantially more than the er- 
rors assigned, we have remeasured them using 
a different technique in order to obtain a new and 
independent set of values. The high efficiency 
for conversion of gamma rays from 7° mesons 
into electron pairs and the ease of recognition 
of electrons in a xenon bubble chamber* have 
been used to separate the various K*-decay 
modes in our experiment. 

The results presented in this Letter are based 
on the analysis of a sample of 6300 K* mesons 
obtained from an exposure of the xenon chamber 
to a separated 700-Mev/c K* beam at the Beva- 
tron. A moderator was used to reduce the beam 
momentum to about 400 Mev/c at the entry of the 
chamber, so that the K* mesons came to rest 
and decayed near the center of the chamber. The 
chamber has no magnetic field. 

The scanning and measurement of the photo- 


graphs led to the classification of the events into 
various categories. The requirements that char- 
acterize each of these categories are given below, 
followed in parentheses by a list of the K*-decay 
modes that are the main constituents of the events 
classified therein: (a) events with more than one 
charged secondary (rT); (b) events with a single- 
charged secondary and three or four electron 
pairs (r’); (c) events with two, one, or zero elec- 
tron pairs and a charged secondary clearly iden- 
tified as an electron by its large radiative energy 
loss or its production of a cascade shower (Ko3). 

The remaining events were classified as fol- 
lows: (d) events with a single-charged secondary 
and two electron pairs identified to be Kz by the 
constraints imposed in the application of energy 
and momentum conservation to K* and 7° decays 
(Ky2); (e) events not in Category (d) which have 
a single-charged secondary and two electron pairs 
(Ky.3,7’); (f) events with a single-charged second- 
ary and one electron pair (K72, Ky3); (g) events 
with a single-charged secondary and no electron 
pairs (K 2): 

The » secondary from a 7* decay at rest has a 
range of only 1.3 mm and hence cannot always be 
reliably identified. Thus no separation was made 


Table I, K* branching ratios, in percent. 








AB for 2% 

change in y 

o b ‘ Present scanning 

Decay mode Birge et al. Alexander et al. Taylor et al. expt. efficiency 
Kote +e 58.5 +3.0 56.9 +2.6 64.2+1.3 0.2 
K oat +n! 27.7+2.7 23.2+2.2 18.6 +0.9 0.1 

TT 

Kye tas y 2.8+1.0 5.9+1.3 2.8+0.4 4.8+0.6 0.3 
K,7e +04 3.2+1.3 5.1+1.3 5.0+0.5 0.0 
t — onten™ 5.6 +0.4 6.8+0.4 5.2+0.3 5.7+0.3 0.0 
tT '— t+ 2n° 2.1+0.5 2.2+0.4 1.5+0.2 1.7+0.2 0.1 








a 
See reference 1. 
See reference 2. 
See reference 3. 
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between 7 and yp secondaries on the basis of ob- 
served 7-u. decays. Since most K* mesons stop 


near the center of the chamber, the 24-cm range 


of the 7* from the Ky2* decay mode is too long in 
a 30-cm diameter chamber to be useful as a meth- 
od of separation of K,2 from other modes. 

We now illustrate in a slightly oversimplified 
way the method of calculating the branching ratios 
from the data. We let N,, N,,+++N, denote the 
numbers of events in the above seven categories, 





B(t)=N,/N, 
= N, 
~ N[C(r’, 4) +C(7’, 3)’ 


B{K 4) =N,/CN, 





B(t’ 





1 
BKK 3) “OR 52) 3) [N./N -C(r’, 2)B(T ) -B(K,,)(1 - 
B(K .)= ! 
12) “CK, 2) +O 


and let N =));N;. We also define C(K;, j) as the 
probability that in the decay mode K;, 7 gamma 
rays are converted and observed in the appropri- 
ate fiducial volume of the bubble chamber, e.g., 
C(K y3> 1) is the probability of observing one gam- 
ma ray from a Ky3. Finally we let Ce be the 
probability of recognition of an electron, by vir- 
tue of its large energy loss or production of a 
shower in the chamber. In terms of the above 
quantities, the branching ratios can be calculated 
from the following relations: 


CICK 42)) 


Ay Og +N) - Crt, IUBEr) CUE, gy BU, 9) = (1 CICK, 9, NB, 3} 


BKK 4) =N,/N -C(r’, 0)B(r’) -CE Ly O)B(K | 5) -C{K 4, 0)B(K .) -(1 -C CK 3 0)B(K 4): 


The determination of the efficiencies C(K;j, j) 
consists of two parts: computation of the conver- 
sion probabilities in the chamber, and estimation 
of the fraction of electron pairs missed because 
of scanning inefficiency. The first part was done 
by a Monte Carlo technique in which the observed 
spatial distribution of K*-decay points and the 
theoretical gamma-ray conversion cross section 
were used to compute the probabilities of forming 
various numbers of electron pairs within the ap- 
propriate fiducial volume. These efficiencies are 
insensitive to assumptions about 7° momentum dis- 
tributions in the three-body decay modes and are 
approximately represented by multinomial distri- 
butions with conversion probability of 80% for a 
single gamma ray. From the Monte Carlo calcu- 
lation, we find that the fraction of K72 events in 
which neither gamma ray from the 7° converts is 
onl¥ 5%. 

The effects of scanning inefficiency were mini- 
mized in the following ways: (a) All pictures were 
scanned twice by technicians or physicists. (b) 
About 30% of the pictures were scanned independ- 
ently by two physicists. (c) Only pictures of the 
best quality, without too many K* mesons, were 
used in the analysis. (d) All K*-decay points were 
required to fall within an appropriate fiducial vol- 
ume appreciably smaller than the bubble chamber. 





An important check on the efficiency was afford- 
ed by a comparison of the experimentally found 
relative proportion of events with one and with two 
electron pairs and those obtained by the Monte 
Carlo computation. From this comparison, we 
conclude that the average scanning efficiency for 
electron pairs is (85+ 2)%. Much of this loss can 
be ascribed to steeply dipping and low-energy 
pairs. An error of 2% in the efficiency estimate 
would change the branching ratios by the amounts 
shown in the last column of Table I. 

The electron recognition efficiency Cg was ob- 
tained by measuring, for a sample of identified 
Ke3 events, both the minimum amount of track 
required to identify the electron and its potential 
path length in the chamber. From these data, the 
electron recognition probability as a function of 
potential path was evaluated and, averaged over 
the potential path distribution, yielded Cz 
= (89+ 2.5) %. 

Small corrections were applied to the results 
obtained as described above for the following ef- 
fects: (a) the presence of a small contamination 
of decays in flight, for which the kinematical Ky2 
tests discussed previously are incorrect; (b) the 
fact that, because of measurement errors, a 
small number of K,,3 events pass the Ky kine- 
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matical tests. 

The final results are given in Table I, which 
for comparison also contains the results of emul- 
sion determinations by Birge et al., Alexander et 
al., and Taylor et al. 

No clear example of a decay mode other than the 
six considered above was found, but it must be 
emphasized that no systematic search aimed spe- 
cifically at the identification of rare modes has 
yet been completed. 

We wish to express our appreciation to Dr. J. D. 
Van Putten and Dr. H. C. Bryant, and to the staff 
of the Bevatron for their invaluable help in the 
course of this experiment. 
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Atomic Energy Commission. 
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INTERACTIONS OF LAMBDAS WITH PROTONS 


Gideon Alexander ,f Jared A. Anderson, Frank S. Crawford, Jr., William Laskar, 
and Lester J. Lloyd 
Lawrence Radiation Laboratory, University of California, Berkeley, California 
(Received September 29, 1961) 


In the course of an experiment to study the as- 
sociated production processes, 


7” +p+A+K° (1) 
and 
17 +p+D°+K®, D°+Aty, (2) 


we have examined subsequent A interactions in 
hydrogen of the types 


A+p+Atp, (3) 
A+p+Z°+p, (4) 
Atp+ dten. (5) 


In this Letter we present results for reactions 

(3) and (4), based on a sample of 5000 A’s pro- 
duced in the Alvarez 72-in. hydrogen bubble cham- 
ber. [We postpone discussion of reaction (5) to 

a later publication.] In addition, we set an upper 
limit to the pion-producing reactions, 


A+p+A+p+7° 
and 
A+p>A+n+n", (6) 


The A “beam” momentum extends from 400 to 
1000 Mev/c, and corresponds to the associated 
production processes (1) and (2) for incident 7~’s 
of 1035 Mev/c. The distribution of A momentum 
versus path length is shown in Fig. 1. It includes 
two cutoff corrections: (a) Single V° decays of 
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K° or A (without subsequent scatter) are rejected 
if the neutral particle travels less than 0.5 cm 
before decaying; (b) only those A interactions 
that occur within 10 cm of the production point 
are accepted. Cutoff (a) serves to exclude spu- 
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FIG. 1. Distribution of A path length versus mo- 


mentum. The arrows indicate the thresholds for re- 
actions (4) and (6). : 
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rious V°’s due to two-prong events. Cutoff (b), 
while including most of the A path length before 
decay, serves to minimize the geometrical es- 
cape correction, and helps to exclude random 
background recoil protons from candidacy. 

Lambda interactions that are followed by the 
charged decay A+p+m~ are detectable whether 
or not the associated K° undergoes charged de- 
cay. Interactions followed by the neutral decay 
A+n+m7° are also detectable and are used, pro- 
vided that the K° decays via the charged mode. 
The kinematics and measuring errors are such 
that there is no difficulty in distinguishing be- 
tween reactions (3) and (4). 

Some of the interactions were found by scan- 
ners in the first scan. However, the main search 
for interactions consisted of (a) a systematic ex- 
amination on the scanning table of all single and 
double V° events which failed to satisfy the kine- 
matics for production via reactions (1) and (2), 
and (b) an examination on the scanning table of 
all single V° events where the decaying neutral 
is identified as a K®. Here one looks for a re- 
coil proton starting from the line of flight of the 
A. 

Recoil protons of less than 0.5 cm are difficult 
to detect, and have been excluded. Our elastic- 
scattering cross section [Reaction (3)] therefore 
does not include small-angle scatters of the A. 
The excluded angular region depends on the mo- 
mentum of the incident A, and is indicated by the 
shaded region in Fig. 2. 

The correction to the integrated cross section 
depends on the shape of the angular distribution. 
The elastic-scattering cross sections include a 
cutoff correction for which uniformity in cos@, 
(c.m.) and Pa, (lab) is assumed. From Fig. 2 
we see that, within the limited statistics, these 
assumptions are satisfied. The correction factor 
amounts to 1.074 for P, (lab) from 400 to 638 
Mev/c, and 1.028 for P, (lab) above 638 Mev/c, 
the threshold for 5° - p production. 

In the cases of the endothermic processes (4), 
(5), and (6), the proton range exceeds 0.5 cm 
for all production angles over the entire range 
of P, (lab), so that there is no corresponding 
cutoff correction. 

Details of the individual events are found in 
Table I. Figure 2 shows the laboratory momen- 
tum of the incident A and the c.m. angle of the 
final hyperon (relative to the incident A) for each 
event. Table II gives our cross-section results. 
The elastic-scattering cross section is given for 
the entire momentum interval, and separately 


for the intervals below and above 638 Mev/c, the 
threshold for reaction (4). The upper limit given 
for the pion-producing reaction (6) is the cross 
section that would correspond to finding a single 
event instead of the number (zero) actually found. 
Our elastic-scattering cross section can be 
compared with the previous result of Crawford 
et al.,’ obtained in substantially the same momen- 
tum interval. Based on four events, they found 
a cross section of 40+ 20 mb, which is not in dis- 
agreement with our result of 22.3+ 5.9 mb. 
Information from study of hyperfragments has 
been used by Kovacs and Lichtenberg? to calcu- 
late the elastic-scattering cross section at 416 
Mev/c (75 Mev) and at 598 Mev/c (150 Mev). 
They chose a phenomenological central potential 
with a hard core and considered only angular mo- 
menta of <2. Their calculated results are 26 and 
21 mb at 416 and 598 Mev/c, respectively, when 
no spin-orbit term is included, and 34 and 32 mb 
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FIG. 2. Details of the interactions. Solid circles 
indicate elastic scattering events [reaction (3)]. Open 
circles represent 2°-p production [reaction (4)]. The 
shaded region corresponds to elastic scattering in which 
the recoil proton has a range of <5 mm. 
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Table I. 


—$$—$—_—_ 


Details of the interactions. 





Initial lab momentum 
of A 
(Mev/c) 


Center-of-mass momentum cos6y(c.m.,) 


Initial Final 





A. At+p—~A+p 


602 103 
722 504 
739 246 
865 515 
818 550 
586 135 
537 508 
834 420 
864 246 
814 534 
602 317 
857 282 
539 021 
724478 


B. A+p—Z°+p 


861205 
586 381 
692 127 





when spin-orbit interaction is included. Our 
measured value of 24.7+ 9.3 mb for the interval 
400 to 638 Mev/c is compatible with either cal- 
culated result. Thus we agree with the predic- 
tion from hyperfragment data, but cannot distin- 
guish between the presence or absence of spin- 
orbit terms. 

In the doublet approximation (DA),’ the >-A 
mass difference is neglected and 5t,0 and A are 
combined into the two doublets, 


N,=(Z*, (A-2°)/v2] and N,=[(A+D°)/V2, D7). 


Both N, and N, have the same coupling to the 
pion field. This pion coupling need not be the 
same as that of the nucleon doublet N,=([p, n]. 
In addition, a new quantum number is introduced 
which prevents mixing of the doublets.* (For in- 


stance, 5~+p+2*+n+127 is forbidden in DA.) 
The mass difference between N, or N, and the 
nucleon N, is not neglected. 

One can then show that 


2[o(Ap + Ap) + o(Ap + D°%)] 


=0(D*p +=D*p)+0(5p +d). (7) 


Of course, Eq. (7) can only be valid at energies 
high enough that the actual >-A mass difference 
is negligible compared with the c.m. kinetic en- 
ergy. For instance, Eq. (7) is certainly inap- 
plicable below threshold for Ap+2%. As a first 
approximation to a test of Eq. (7), we compare 
our results for the reactions on the left side of 
Eq. (7) to those of Stannard‘ for those on the right 
side. We consider the Ap + Ap cross section only 


Table II, Summary of cross-section results. 





Reaction P , (lab) 


(Mev/c) 


No. of 
events 


A path length 
(cm) 


Cross section 
(mb) 





A+p—At+p 
A+p—A+p 
A+p—Atp 
A+p— Z°+p 


A+p+n° 
ang oe 


400-1000 
400- 638 
638-1000 
638-1000 


880-1000 


1 


18770 

8700 
10070 
10070 


1986 


22.3+5.9 
24.7+9.3 
20.4+7.7 

8.5+4.9 


<14+14 
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for Py, (lab) above 638 Mev/c, the threshold for 
reaction (4). Our range of total c.m. energy cor- 
responds to a >-p c.m. momentum range of 0 (at 
threshold) to 314 Mev/c [for P ‘4 (lab) of 1000}. 
Stannard’s results correspond to a >-p c.m. mo- 
mentum range of 0 to about 600 Mev/c. 

We find 


2[o(Ap + Ap) + o(Ap + 5% ]=57+18 mb. 
Stannard finds 
o(*p + =*p)+ o(S pp +d p)=48+17 mb. 


Within the statistical errors, the prediction [Eq. 
(7)] of the doublet approximation is well satisfied. 
In the theory of global symmetry (GS),° the as- 
sumption (in addition to those of DA) is made that 
the pion coupling of N, and N, is the same as that 

of the nucleon N,, and the mass difference be- 
tween N, or N, and N, is neglected. One can then 
predict N,-N, or N,-N, interactions from N,-N, 
(nucleon-nucleon) interactions. de Swart and 
Dullemond have used potentials that fit avail- 
able nucleon-nucleon scattering data, in linear 
combinations prescribed by GS, to predict cross 
sections for hyperon-nucleon interactions.® They 
take exactly into account the >-A mass differ- 
ence in the kinematics (but neglect the © -nucleon 
mass difference). Their predictions for A-p 
elastic scattering and for 5°-p production [re- 
action (4)] are shown in Fig. 3, together with our 
experimental results. Within statistics, our 
cross sections agree with the predictions of GS. 

de Swart and Dullemond also predict angular 
distributions for A-p elastic scattering at 475, 
640, and 815 Mev/c. These momenta cover our 
range of P, fairly well (see Fig. 1). With our 
limited statistics we can check at most a single 
parameter of the angular distribution. We ex- 
amine our ratios for forward-to-backward scat- 
tering (F/B) and for polar -to-equatorial scatter - 
ing (P/E). Since we see no strong momentum 
dependence in the elastic-scattering cross sec- 
tion, and since the predicted angular distributions 
do not vary drastically over our momentum range, 
we average with equal weight the theoretical pre- 
dictions for F/B and P/E and compare them with 
our entire sample of 14 events. 
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FIG. 3. Comparison of our experimental results 
with the predictions (based on global symmetry) of 
de Swart and Dullemond.* The upper curve and the 
solid circles represent elastic scattering. The lower 
curve and open circle correspond to r°-p production. 


For F/B we find 7/7=1.0+0.5. The prediction 
of GS is 1.5. 

For P/E we find 9/5=1.8+1.0. The GS predic- 
tion is 2.0. 

It is a pleasure to acknowledge the guidance 
and encouragement of Professor Luis W. Alvarez. 
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1™*+p AND 1” +p TOTAL CROSS SECTIONS FROM 8 TO 20 Bev/c* 


S. J. Lindenbaum, W. A. Love, J. A. Niederer, S. Ozaki, J. J. Russell, and L. C. L. Yuan 


Brookhaven National Laboratory, Upton, New York 
(Received October 12, 1961) 


Previous experiments’’? by von Dardel et al. 
have reported measurements of the 7*++ p and 7~ 
+ p total cross sections over the incident pion mo- 
mentum range of 4-10 Bev/c. These investigators 
have found that both the 7* +p and 1” +p total cross 
sections decrease slowly and monotonically (by 
=10%) going from 4.5 Bev/c to 10 Bev/c incident 
pion momentum. 

At 10 Bev/c, o;(7~ +p) =27.0+ 0.4 mb, o¢(1* +p) 
=25.2+0.4 mb, and o¢(7~ +p) - of(7* +p) =1.84 0.6 
mb, where the subscript ¢ denotes total cross sec- 
tion. As those authors concluded, the data do not 
satisfy the necessary conditions for application 
of the Pomeranchuk theorem*® which is that both 
oz(n*+ +p) and o;(1~ +p) become and remain con- 
stant with increasing energy. As a matter of fact, 
neither cross section appears to be consistent 
with this condition. 

In a recent publication* the authors of this com- 
munication have investigated the antiproton-proton 
and the proton-proton total cross section from 4 
to 20 Bev/c. 0,(p+p) appeared to become and re- 
main constant within the errors from 10 to 20 Bev/c 
[o,(p+p)~39.5+1 mb], but the total antiproton-pro- 
ton cross section continued to decrease monotoni- 
cally from 10 Bev/c where o,(b +p)= 58+ 4 mb to 
20.3 Bev/c where o;()+p)=48+4 mb, and hence 
the Pomeranchuk theorem could not be verified 
by this investigation. Therefore it appeared to be 
of considerable interest to extend the investigation 
of the o,(7~ +p) and the o;(7* +p) to the so far un- 
investigated energy range of 10-20 Bev/c. The 
1~-n+ particle-antiparticle pair is especially in- 
teresting to compare since, in addition to the gen- 
eral Pomeranchuk theorem,’ there is another spe- 
cific argument which can be applied to the beha- 
vior of the difference o,(7~ +p) - o;(* +p) at high 
enough energies which cannot be applied to the be- 
havior of the difference o,() +p) - 0;(p +p) or the 
difference o,(K +p) - o,(K*+p). This argument, 
which was proposed by Okun and Pomeranchuk,® 
can be summarized as follows. At high enough 
energies one would expect the charge-exchange 
cross section to become small compared to the 
total inelastic cross section, since charge ex- 
change represents only one inelastic channel of 
what one might expect to be an arbitrarily large 
number of possible inelastic channels at sufficient- 
ly high energy. Therefore the T = 3 state complex 
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scattering amplitude must approximately equal 
the T = 4 state complex scattering amplitude, since 
the charge-exchange amplitude is proportional to 
the difference of these two amplitudes. This then 
implies that the 7” +p total elastic scattering 
cross section equals the 7* +>) total elastic scat- 
tering cross sections and hence also that 


o,(a~ +p)= o(a* +p). 


Of course the above argument does not depend 
on the constancy of cross sections at high energy 
which is a necessary requirement for the Pom- 
eranchuk theorem, and also does not depend on 
the validity of the forward dispersion relations 
which are employed to prove the theorem. 

Weinberg’ has recently presented a generalized 
Pomeranchuk theorem which relaxes the condi- 
tions required to prove that o,-0. +0 as E +~; 
o, is the particle total cross section and o- is 
the antiparticle total cross section. The essen- 
tial assumptions are: (a) that the difference 08 
-o_(E) does not change sign an infinite number of 
times; (b) that the cross sections o, do not grow 
faster than (InE)’* as E+, Then, using the for- 
ward dispersion relations, he demonstrates that 
o,-0_.70 as E +~. 

The present experiment to determine o,(a* + pb) 
and o,("~ +p) from 8 to 20 Bev/c was run during 
the same period and with the same experimental 
arrangement as the previously reported work on 
o,(P +p) and o,(p +p),* the only difference being 
that the gaseous Cerenkov counter was tuned to 
a pressure which selected pions rather than pro- 
tons, or antiprotons. The previous method of us- 
ing six counters of different area to extrapolate 
the measured transmission to zero area was 
again employed. The measured transmissions 
and hence also the values of the cross sections 
were corrected for the following effects: 

(a) Due to the m-p and K-y, decays there is a mu- 
on contamination of the beams. Using the known 
pion and K-meson production spectra and angular 
distributions from previous work,® the known K- 
meson and pion lifetimes, and the beam setup 
geometry and magnetic geometry, we were able 
to estimate the muon contamination of the pion 
beam to range from 2 to 3%, with an error of 
+1%; this correction was applied to the data. The 
uncertainty in the difference of oO, +p- o,* +9 
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FIG, 1. 7 +p and7 +p total 
cross sections. The solid lines 
are Fit 2’ and the dashed lines 
are Fit 2 (see text). The mo- 
mentum of the system is defined 
to ~+2.5%, and the absolute scale 
of the momentum is uncertain by 
2h. 
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due to muon contamination is only a fraction of a 
percent. 
(b) The electron contamination of the beam was 


estimated to be ~0.3 % and a correction was made. 


(c) Accidentals mainly in the end counter array 
were monitored during the run and the appropri- 
ate correction (~1%) to the cross sections was 
made with an uncertainty of a fraction of 1%. 

(d) The proton contamination correction to the 
cross section is less than 1% below 15 Bev/c and 
rises to approximately 3.5% at 20 Bev/c. The 
uncertainty in the cross section due to this cor- 
rection is less than 1%. 

The results are shown in Fig. 1. The errors 
on each point predominantly represent point-to- 
point errors. In addition we estimate that the ab- 
solute scale uncertainty may be as much as 2% 
(i.e., 0.5 mb); however, this absolute uncertainty 
has a negligible effect on the difference o;(1~ +p) 
~og(1* +p). 

From Fig. 1 it is clear that our results agree 
very well with those of von Dardel et al. at the 


8 10 l2 14 
MOMENTUM P (Bev /c) 


lower momenta. Considering the errors, it ap- 
pears that there is no convincing evidence from 
the present data for other than a smooth slow vari- 
ation of the total cross section with energy. There- 
fore we have attempted three least-squares fits 

of the momentum dependence of the total cross 
sections from 8 to 20 Bev/c. Fit 1: o,=a; Fit 2: 
o,=a+b/p; Fit 3: o,=a+b/lnp; where both a and 
b are constants. It was found that Fit 2 has the 
lowest x”, and on inspection of the lower energy 
data of von Dardel et al. it appears that Fit 2 and 
Fit 3 are the only fits capable of reasonably fit- 
ting the data from 4.5 to 20 Bev/c. 

Fit 2’ is the Fit 2 momentum dependence o;=a 
+b/p but applied to the sum of our data and the 
data of von Dardel et al. (i.e., a fit from 4.5 to 
20 Bev/c). Both Fit 2 and Fit 2’ are shown in 
Fig. 1. It is clear that this functional dependence 
can be made to fit both our own and the lower en- 
ergy data reasonably well. The parameters a and 
b for these fits are given in Table I. 

The calculated differences A=0;@~ +) -o;@*+p 


Table I. Fit of o,=a+b/p, where 0; is in mb and p is pion lab momentum in Bev/c. 





Fit 2 
ba 


Fit 2’ 





22.59 +0.63 
25.26 +0.46 


+ 
O(n +p) 
O(n +p) 


19.39+8,.21 
13.03 5,27 


22.26 0.33 
24.37+0.29 


25.10 +2.83 
24.94 42.65 
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for all of the fits 1, 2, 2’, and 3 show that A is 
probably constant or at best slowly decreasing 
from 4.5 to 20 Bev/c. On the other hand, the to- 
tal cross sections decrease monotonically by 
about 10% from 4.5 Bev/c to 10 Bev/c, and the 
additional decrease from 10 Bev/c to 20 Bev/c 

is only about 4%. If the assumed functional de- 
pendence o=a+b/p were to be correct for higher 
energies, our results imply that the limiting val- 
ues of both total cross sections at very high ener- 
gies would decrease from the values at 20 Bev/c 
by about 5% but the same order of difference (A 
~2 mb) would still remain. We also found this to 
be the case for the other fits. 

Although extrapolation of the observed momen- 
tum dependence of the cross section in the ob- 
served energy range to much higher energies is 
not to be taken very seriously, it is, on the other 
hand, strongly implied by our data that at 20 Bev/ 
c we are still very far from the energy region 
where the Pomeranchuk or Weinberg theorems, 
etc., can be applied or even apparently verified. 
Furthermore the substantial difference of (1 ~ +p) 
- of(a* +p) and its slow decrease with energy that 
we have observed also imply that the number of 
inelastic channels which are effective is much 
too small and is growing much too slowly with 
energy to allow us to believe that we are ap- 
proaching the necessary energy to demonstrate 
the high-energy limit theorems considered. 

In view of our previous work on antiproton- 
proton and proton-proton total cross sections’ 
and other previous work on K*-proton and K~- 
proton total cross sections, we now believe that 


—$—__ 


it is likely that even a possible experimental ver- 
ification of the high-energy limit theorems (i.e., 
o,-0_+0 as E ~~) will have to await the next 
large increase in proton accelerator energy. 

The authors wish to acknowledge the coopera- 
tion of the AGS staff and the operating crew dur- 
ing these experiments. The authors also wish to 
acknowledge the assistance of J. Elias and C. Sin- 
clair during these experiments. 
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EXISTENCE OF TWO T =0 VECTOR MESONS* 
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We wish to interpret the various recently dis- 
covered vector meson states within the frame- 
work of the vector theory of strong interactions’? 
(VTSI) proposed nearly two years ago. In partic- 
ular, we concentrate our attention on the two dis- 
tinct T=0 vector mesons unambiguously predicted 
by VTSI. The following are some of the major 
points we would like to make: 

(1) The 550-Mev peak in the 7*1~7° mass plot 
for the reaction 


tt+dont+n +97°+2p 


studied by the Johns Hopkins— Northwestern group* 


may well be due to a second T =0 vector meson 
state (tentatively called n°) with exactly the same 
quantum numbers as the experimentally well es- 
tablished w° at 785 Mev. 

(2) In the language of VTSI, the w® is to be iden- 
tified with the vector meson coupled to the baryon- 
ic current whereas the 7° with a lower mass is to 
be identified with the vector meson coupled to the 
hypercharge current. 

(3) Most of the hard-core effect between two 
nucleons should come from the exchange of an w® 
rather than from the exchange of an 7°. 

(4) The conjectured n with a mass lower than 
the p (T=1, J=1- two-pion resonance) mass is 
extremely helpful in understanding the nucleon 
structure. 

(5) Generalized Clementel-Villi-Fubini type 
analyses applied to the observed nucleon form 
factors suggest that 120% of the isoscalar charge 
is due to n, and that 120% of the isovector charge 
is due to p; hence, using the argument of Gell- 
Mann and Zachariasen,* we infer that the univer- 
sality relations implied by the couplings of the 
vector mesons to the appropriate conserved cur- 
rents may hold approximately even at s=m p” 

My» 

Some time ago a theory of strong interactions 
was proposed which unambiguously predicts the 
existence of one J =1 vector meson (p-wave two- 
pion resonance) and two T =0 vector mesons (T 
=0, J=1- three-pion resonances). Recent mN and 
bp experiments®*® have conclusively established 
that there indeed exist a T=1 vector meson p 
with a mass of 750 Mev and a width of ~80 Mev 
and a T=0 vector meson w° with a mass of 785 


Mev and a width <30 Mev. Moreover, Pevsner 
and collaborators® report that, in the reaction 


tt+dont+n~ +17°+2p, (1) 


in addition to the clearly identified w° peak (whose 
mass agrees well with that of Magli¢ et al.°), there 
is an indication for another peak in the Q-value 
plot for 7*1~7° which might be associated with a 
three-pion resonance with a mass of ~550 Mev 
and a width <30 Mev. The statistics are still 
meager, and the effect may disappear later.’ 
Neither the isospin nor the spin-parity assign- 
ment has been made. Yet in the present note we 
take this state seriously, and call it n°. More- 
over, we assume that 7° is a T=0 vector meson 
that decays strongly into three pions with T =0, 
J=1-, just like w®. (As the 7 mass is not far 
above the three-pion threshold, the two-body 
electromagnetic decay mode 7°+y might also be 
detectable; if we use the method of Gell-Mann 
and Zachariasen* to compute the decay rate for 

n° +7°+y from the 7° lifetime, we obtain about 
0.03 Mev for the partial width provided that fy /4n 
~2.) 

Let us recall that the two T =0 vector mesons 
of VTSI are coupled linearly to the two exactly 
conserved T =0 currents of the strong interac- 
tions—the hypercharge current and the baryonic 
current. The following question naturally arises: 
Is w°(n°) coupled to the baryonic current or to the 
hypercharge current? Note that here we have two 
particles or resonant states that have identical 
isospins and identical spin parities; moreover, 
their widths seem narrow in both cases. The only 
difference is that the w mass is 230 Mev higher 
than the 7 mass. Yet according to VTSI, their 
roles in strong-interaction physics are quite dis- 
tinct. For instance, if the pseudoscalar mesons, 
nm and K, are to emerge as tightly bound states of 
NN and NA systems, glued by a heavy neutral vec- 
tor meson as suggested by Teller and others,’* 
the coupling of the baryonic vector meson (BB) 
of reference 1) must be much stronger than the 
coupling of the hypercharge vector meson (BY ) 
of reference 1); otherwise there would be a very 
tightly bound state of a K and an N. 

In principle, it is possible to settle this ques- 
tion by studying the analytic structure of the KN 
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scattering amplitudes as functions of momentum 
transfer. For instance, if n° (but not w®°) is the 
vector meson coupled to the hypercharge current, 
both the T=1 KN amplitude and the T=0 KN am- 
plitude must have poles at ¢=m,,” with equal resi- 
dues, but no poles should be present at ¢=m,,’. 
In practice, however, this will not be feasible for 
a long time to come. A more practical way is to 
study the effects of n° and w® on nuclear forces 
at short distances, keeping in mind that the coup- 
ling of the baryonic vector meson to the nucleon 
must be much stronger than the coupling of the 
hypercharge vector meson to the nucleon. Just 
as the baryonic vector meson is responsible for 
most of the attraction necessary to bind a baryon 
and an antibaryon to form a7 or aK, the same 
particle must also be responsible for most of the 
observed short-ranged, very strong repulsion be- 
tween two nucleons. In this connection we should 
recall an important observation made by Breit,° 
who also advocated the neutral-vector-meson ap- 
proach to nuclear forces to account for the repul- 
sive core and the spin-orbit force in a unified man- 
ner.’° He demonstrates that with a vector coupling 
constant of the order of 10 and with a vector meson 
mass ~4m,, the repulsive core radius would be- 
come too large in the sense that the central force 
in the intermediate region (distance ~ 0.75/m,) 
would be completely dominated by the strong re- 
pulsion due to the vector meson. Although a rig- 
orously quantitative discussion of Breit’s argu- 
ment must await a more reliable estimate of the 
two-pion exchange contribution to nuclear forces 
than is now available, there is no doubt that be- 
tween w with m,)=5.9m_ and 7 with mp=4.0mz,, 
the w is the more likely candidate for the T =0 
vector meson coupled strongly to the baryonic 
current. In the language of reference 1, this 
means that w° is the vector meson associated 
with the B, ‘B) field while n° is the one associated 
with the By, ") field.™ 

It is encouraging that, in both pp collisions and 
m*d collisions, w mesons show up much more con- 
spicuously than 7 mesons. In reaction (1) at p,{tab) 
=1.23 Bev/c, the cross section for w° production 
seems to be about four times as large as that for 
n° production.? In p+p+5z, 77, there is hardly 
any evidence for n’s.° Thus, contrary to Chew and 
Frautschi,’* all strong interactions do not seem as 
strong as possible; rather some strong couplings 
such as the coupling of w° to NN seem much strong- 
er than other strong couplings such as the coupling 
of n° to NN, as conjectured earlier.’ 

We now examine the nucleon structure problem 
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in the light of the existence of p, 7, and w. Elec- 
tron-proton scattering experiments carried out at 
Stanford and Cornell** have revealed that the neu- 
tron charge cloud has a fairly large positively 
charged “fringe.” In a more sophisticated lan- 
guage, the average mass state responsible for the 
isoscalar form factor must be lower than that re- 
sponsible for the isovector form factor. Applying 
Clementel-Villi* type fits to the nucleon form fac- 
tors, it was concluded by Bergia et al.’® (BSFV) 
that, if a simple resonance picture (which has its 
origin in the pioneer work of Nambu"*) holds, the 
T=1, J=1- two-pion resonance proposed by Fra- 
zer and Fulco’”’ must have a higher mass than the 
T =0, J=1- three-pion resonance proposed by 
Nambu?® and Chew.”® Thus the simple resonance 
picture of BSFV would run into trouble if the w 
which has a higher mass than the p were the only 
T =0 vector meson (or the only 7 =0, J=1- three- 
pion resonance). 

It is evident that the introduction of a second neu- 
tral vector meson with a mass lower than the p 
mass will help remedy this situation. Yet one may 
still argue that if n° is coupled less strongly than 
w° to NN, then the 7 meson effect might not help 
too much. This, however, is not necessarily the 
case. It can be shown, using an argument origi- 
nally given by Gell-Mann and Zachariasen,* that 
if the “universal” coupling constant of 7° to the 
hypercharge current defined at s=0 does not differ 
appreciably from the nNN coupling constant defined 
at s=m,?, then there must be a substantial contri- 
bution to the isoscalar charge form factor from the 
one 7 state, the crux of the matter here being that 
both 7° and the “isoscalar photon” are coupled to 
the same conserved hypercharge current. More- 
over, in a theory which is sufficiently symmetric 
between N and =, it is easy to show that the transi- 
tion “isoscalar photon” + w° must be forbidden to 
the extent that the N= mass difference could be ig- 
nored. (Recall that the baryonic current is even 
under N= =, while the hypercharge current is 
odd.) 

But let us work backwards, so to speak, and 
look at the experimental data first. If there are 
two T =0 vector meson states, it is natural to 


modify the BSFV formula in the following way: 
2 2 
w 

5+(1 “4, -a ). (2) 
w 


am am 


s ae w 
Fy ST 





Crudely speaking, Gp, @,, and (1-ay - Ay) rep- 
resent the fractions of the isoscalar charge contr! 
buted by the one 7 state, the one w state, and all 
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other (hopefully very massive) states, respectively. 
if we regard m,”=16m,” and m,,? =32 m,” as ex- 
perimentally determined quantities, the formula 
contains only two independent adjustable param- 
eters, one of which is more or less determined 

by the requirement (71*),”* =0.8 x10" cm, (71”)/” 
=0. Crude numerical estimates show that the ob- 
served data can be adequately reproduced by’® 


a =1.2, a =-0.7. 
n w 


On the other hand, for F,”(q?) the original BSFV 
form suffices so ne we wee 


a= 


ee a . 
jar *) (3) 


just as before; the data require 


a =1.1-1.2, m_?=19m ?-22m ?. 

p p 7 7 
(The slight “discrepancy” between this value for 
m,* and the “observed” value mp” = 29m,” is some- 
what puzzling.) 

Gell-Mann and Zachariasen* have noted that if p 
and n are coupled to the conserved isospin and the 
hypercharge current, respectively, and if the bare 
masses of p and 7 are infinite (or very large), 
then the constants @p and Gy are directly related 
to the form factors of the vector type couplings of 
the vector mesons to the nucleon at zero momen- 
tum transfer, or, equivalently, to the ratios of 
the “universal” coupling constant fp /4n or fn /4n 
at zero momentum transfer to the _ coup- 
ling constant measurable at s = =Mp° , my? denoted 


by fonN* /4n, Snnn’/47: 

a= U/F oy w's)],. 07 Sonn’ , 
(4) 

«, = [L/F py wis) - 0 “fanny 

It is extremely gratifying that the observed values 

of a» and @p are so close to unity. We do not yet 

understand the deep reason for this, but in view of 

(4), relations such as @p)~1 and a,~1 seem much 

more plausible in a theory in which the vector mes- 

on states coupled to the conserved currents are 

introduced in the very beginning, rather than in a 

theory in which these rather narrow resonant 

States emerge dynamically, if they do at all, ina 

mysterious “self-generating” manner 4 la Chew.”° 

Should the a@’s for other charged particles also 

turn out to be approximately equal to unity, then, 

ven within the framework of dispersion theory 

in which all energy momenta are on the mass 

Shell, it would become meaningful to talk about 


the universal coupling of p(n) to the particles bear- 

ing isospins (hypercharges). Initial steps along 

this line have already been made in comparing 

the effect of p on ™N scattering (proportional to 

pane an) With the width of p (proportional to 
as). Mya 22 It appears likely that the notions of 

cnuegnaity and conserved vector currents are 

important elements in our quantitative, as well 

as qualitative,’ understanding of the dynamics of 

strong interactions. 
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EXPERIMENTAL TESTS OF MACH’S PRINCIPLE* 


R. H. Dicke 
Palmer Physical Laboratory, Princeton University, Princeton, New Jersey 
(Received September 11, 1961) 


In this note it will be shown that, contrary to the 


suggestion of Cocconi and Salpeter,’ the extremely 
precise null result of the experiments of Hughes, 
Robinson, and Beltrow-Lopez? and Drever' is to 
be expected, resulting from an application of 
Mach’s principle. 

According to Mach’s principle, as formulated 
by Bishop Berkeley,* Mach,° and Sciama,® the in- 
ertial forces experienced in an accelerated labo- 
ratory are gravitational, having their origin in 
the distant matter of the universe, accelerated 
relative to the laboratory. Because of the tensor 
character of the gravitational-inertial field, it 
should exhibit tensor polarization properties. In 
particular, as suggested by Cocconi and Salpeter, 
because of the flattened rotating mass distribution 
of our galaxy, the inertial reaction having its ori- 
gin in this mass distribution should exhibit some 
spatial anisotropy. This should appear in the for- 
malism as a tensor inertial mass. Cocconi and 
Salpeter suggested that if Mach’s principle were 
valid, the effects of this tensor inertial mass 
would appear as a spatial anisotropy in certain 
experiments. Several experiments, designed to 
test Mach’s principle in this way, have been per- 
formed or analyzed.*"**? By far the most accurate 
has been that of Drever.* 

It will be shown that the experiments do not rep- 
resent a test of Mach’s principle in the manner 
suggested by Cocconi and Salpeter. On the con- 
trary, and in agreement with the requirements 
of Mach’s principle, the experiments show that, 
with great precision, the anisotropy of the inertial 
mass is universal, the same for all particles. 

Expressed relativistically, the suggestion of 
Cocconi and Salpeter* is that the four -momentum 
of a particle can be written as 















P.=m.u? (1) 
= 


where u? is the four-velocity of the particle. In 
the absence of a gravitational field, the particle 
is assumed to obey the equations of motion, 


aP /ds=F , (2) 
1 t 


where F; is the four-force acting on the particle 
and as usual F; satisfies the condition 


Fu’ =0. (3) 


As Mach’s principle associates the inertial reac- 
tion with the matter distribution in the universe, 
an anisotropy in the inertial mass should be uni- 
versal, the same for all particles at the same 
space-time location, for all particles would see 
the same mass distribution. With the assumption 
that the inertial reaction is universal, the tensor 
mj; can be expressed as 


ms, = mh 9 (4) 


where fij is a universal tensor field (dimension- 
less). 

A serious objection to Eqs. (1)-(4) can be raised. 
Equations (2) and (3) are generally not consistent 
with the constraint condition, 

ii 
8, u ,=i. (5) 

Consistent equations can be obtained from a var- 
iational principle. We note that to generate a mo- 
mentum such as Eq. (1), linear in a four-velocity, 
a first condition to be satisfied is that the Lagran- 
gian of the particle should be a function of the in- 
variant 


(d5/ds)* =f,piu! (6) 


Equations of motion of a free particle (gravitation- 
al forces only) are obtained from the variational 
principle, 


0=6fLds, (7) 


where variations of the coordinates are to be tak- 
en subject to the constraint, Eq. (5). With the as- 
sumption that Sij and &jj are not equivalent, the 
resulting expression for the four-momentum is 
linear in some four-velocity only for the unique 
choice, 

i jv2 


L=m(f,wu (8) 


With this choice, the constraint condition, Eq. (5), 
is satisfied by the equations of motion automatical- 
ly, without introducing the condition of constraint 
in the variational calculation. As a result, the 
equations of motion do not contain 8ij explicitly. 
Substituting Eq. (8) in Eq. (7), the equations of 
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motion of a free particle are 


ul m. wea” 
ij 1 jk,i 


ds Fem) ve 


m 
0. (9) 





-_ 2 ~ = 
UF jVe 


These are most conveniently expressed, intro- 
ducing fij a8 a new metric tensor, by defining 


dv’ « f,gexiax’, (10) 


i’ =dx' /d3. (11) 


With these substitutions, Eq. (9) becomes 


—jk 


d j 
)- am, 


(mi 


a "ij (12) 


a 
"2 
The resulting particle trajectories are geodesics 
of the new geometry, with fjj a8 metric tensor. 
The limiting trajectories of particles with infinite 
energy are null geodesics of the new geometry. 

Inasmuch as 8ij does not appear explicitly in the 
classical equations of motion of a particle, the ap- 
propriate quantum mechanical wave equations, 
giving equations of motion of expectation values 
equivalent to the classical equations, are construct- 
ed by employing Sij as the metric tensor. For ex- 
ample, the appropriate Lagrangian density for the 
Klein-Gordon wave function is 


wn? Do Y., ~tm*9*. (13) 
»t 5J 

It may be noted with the assumption made above, 
that the inertial reaction is universal, the same 
for all particles including photons and pions, the 
metric tensor &jj appears nowhere in the formal- 
ism. In fact, for the geometry defined by meas- 
urements in the usual way with real rods and 
clocks, jj is the metric tensor.® 

It should be noted also, that because of the uni- 
versal character of the inertial anisotropy, being 
present in the same way for all particles (or 
fields), the spatial anisotropy is unobservable lo- 
cally. The easiest way to see this is to note that 
a coordinate system can always be chosen in such 
a way as to cause fj; to be locally Minkowskian 
with vanishing first derivatives. For this coordi- 


————_ 


nate system, the anisotropy of inertial mass ap- 
pears explicitly nowhere in the equations. 

While, with these assumptions, inertial anisot- 
ropy is not locally observable, the fact that the 
geometry defined by real rods and clocks is non- 
Euclidean may be interpreted as due in part to 
the effect on rods and clocks of the anisotropy of 
the inertial mass of the elementary particles 
which comprise the rods and clocks. 

It is concluded finally that the extremely accu- 
rate null result of the experiment of Hughes et al. 
does not cast doubt upon the validity of Mach’s _ 
principle. On the contrary, and in accordance 
with the requirements of Mach’s principle, this 
important experiment shows, with great precision, 
that inertial anisotropy effects are universal, the 
same for all particles. 
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8what then is the significance of the redundant tensor 
&;;° As was discussed recently [C. Brans and R. H. 

cke, Phys. Rev. 124, 925 (1961)], the choice of units 
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ERRATA 





K,° DECAYS AND INTERACTIONS. D. Luers, 
I. S. Mittra, W. J. Willis, and S. S. Yamamoto 
[Phys. Rev. Letters 7, 255 (1961)]. 


On page 257, the first line at the top of the page 
reads: “Using the AJ= rule, he predicts a 
= -0.0109+ 0.022. We find a=0.0171+0.0065.” It 
should read: “We find a = -0.0171+ 0.0065.” 


FIRST -ORDER TERRESTRIAL ETHER DRIFT 
USING THE MOSSBAUER RADIATION. Martin 
Ruderfer (Phys. Rev. Letters 5, 191 (1960)]. 


The theory of the proposed experiment is in- 
complete because the first-order variation of 
frequency with motion is not considered. In the 
Fitzgerald-Lorentz contraction theory, the fre- 
quency of a moving clock is f,(1 -w*/c?)”*, where 
f, is the frequency when the clock is at rest in the 
ether and w is its velocity through the ether. For 
the radiator at the center of the turntable, w =v 
and the frequency, f, of the radiator is 


f=fo(l -v?/c?)™. 


For the absorber, w=|¥+@g,| and the frequency, 
fq, of the absorber is 


fe =f, - (ws? +v? - 2wsv sind) /c? . 


The relative frequency shift due to clock motion 
is, to first order, 


\Af/fl a” (f, - f)/f =wsv sind /c?. 


It is readily verified that this is also obtained 
when the radiator is not at the center of the turn- 
table. This frequency shift is equal to the fre- 
quency shift found due to phase shift in rotation. 
The contraction theory therefore predicts a null 
effect as does relativity theory for a one-way 
rotating terrestrial ether drift experiment. The 
proposed experiment is not a crucial experiment 
for deciding between the two theories. To con- 
form to both theories, the proper interpretation 
of the predicted null result is that detection of 

an ether is precluded as required by the special 
theory of relativity and that existence of an ether 
is permitted as required by the contraction theory. 





